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INTRODUCTION 



In monolithic semiconductors with a single bandgap 
the ultimate efficiency, only due to limited absorption 
and to thermalization processes, is around 42% (for the 
black-body radiation at Ts = 6000 K^). The absorption 
of photons is incomplete because photons of energy below 
the bandgap width are not absorbed by the system. On 
the other hand, thermalization causes the surplus energy 
of electrons above the bottom of the conduction band 
or holes below the top of the valence band to dissipate 
in thermal contact with the crystal lattice. In ultimate 
efficiency calculations the utiHzed power is determined by 
the flux of absorbed photons and the width of the energy 
gap. 

The efficiency taking into account recombination pro- 
cesses and the effect of thermodynamic losses is referred 
to as detailed balance efficiency. The recombination rate 
mainly depends on the temperature of the system. At 
ambient temperature only a minor part of excited carri- 
ers undergo radiative recombination to the ground state. 
Thermodynamic losses are a consequence of the fact that 
the chemical potential diflerence (for thermalized conduc- 
tion band electrons and valence band holes) is lesser than 
the gap widtb^. Thus, the external voltage, even with- 
out load, is lower than that resulting from the gap width. 
With progressively increased load the current grows to 
saturate quickly, with a concurrent reduction of voltage 
at cell contacts. In the determination of detailed bal- 
ance efficiency the utilized power is defined as the maxi- 
mum power output at the optimum load. The theoretical 
limit of detailed balance efficiency is around 31% (for the 
black-body radiation at Ts — 6000 K}i. 

The above-mentioned limitation of efficiency is referred 
to as the Shockley-Queisser limit. Many attempts have 



been made recently to find mechanisms and systems that 
would allow to overcome this constraint. Devices in 
which the Shockley-Queisser limit is outperformed are 
referred to as third-generation photovoltaic systems^i^. 

Photon absorption can be increased by using tan- 
dem cells - cascades of solar cells of successively nar- 
rowing bandgapsS*!. A similar effect can be obtained 
with solar-energy converters based on semiconductor 
superlattices^ii^ with multiple bandgap a^^i^^i^^i^^i^^ . 
Multiple-bandgap systems also allow to reduce thermal- 
ization losses. 

Another method of increasing the efficiency is based 
on increasing the photocurrent or the photovoltage. The 
photocurrent can be increased by generation of extra 



carrier* 
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in the inverse Auger effeclji^. Otherwise 
lost in thermaHzation, the energy of highly excited car- 
riers is thus utilized. The effectiveness of this process 
depends on carrier localization and is higher in dot and 
rod superlattices^. The photovoltage can be increased 
in systems operating with hot electrons2ii2^i2^^ in which 
energetically selective electrodes can capture highly ex- 
cited carriers before they attain the thermodynamic equi- 
librium with the lattice through thermalization. 

In this paper we determine the band structure of a 
two-dimensional GaAs/AlGaAs superlattice in the effec- 
tive mass approximation^. Both GaAs and AlAs are 
well-known materials used widely in the fabrication of 
semiconductor heterostructures. Both compounds crys- 
tallize in the zinc blende structure and are character- 
ized by similar lattice constants. In our calculations we 
use the envelope function approximation, which is widely 
applied for the description of conduction band bottom 
and valence band top for zinc blende crystallographic 
structures. We also omit the effect of stress on the ma- 
trix(AlGaAs)/rod(GaAs) interface. AlGaAs heterostruc- 
tures have a relatively wide energy gap (1.4 eV - 2.1 eV) 



which makes them potentially suitable for intermediate 
band solar cells, where a wider gap (e.g. in comparison to 
the bulk Si) is required for achievement of high efficiency 
of solar radiation conversion. 

In a wide range of structural parameter values mini- 
bands in the conduction band as well as those in the va- 
lence band overlap to a large extent. Only the lowest con- 
duction miniband, narrow and quite strongly localized in 
potential wells formed by GaAs rods, is distinctly set 
apart^. Higher minibands, which tend to He above the 
matrix potential, are wide and overlap to form a continu- 
ous block. Even though some isolated minibands can be 
distinguished in valence band wells, the spacing between 
them is very narrow. Therefore, in our estimation of ul- 
timate efficiency the band structure of the solar-energy 
converter is assumed to consist of a continuous valence 
band and a continuous conduction band with a single in- 
termediate band formed by the lowest conduction mini- 
band. 

We determine the ultimate efficiency of a solar-energy 
converter of band structure as defined above, with the 
effects of incomplete photon absorption and carrier ther- 
malization taken into account. The ultimate efficiency 
is examined versus structural parameters of the super- 
lattice; these include the filling fraction, the superlattice 
constant, the rod geometry and the effective potential in 
the matrix material. 

The paper is organized as follows, in the next sec- 
tion we discuss the geometry of the GaAs/AlGaAs super- 
lattice under consideration and present the assumptions 
made in the effective mass approximation for conduction 
band electrons and for interacting heavy and light valence 
band holes. Then we present the plane- wave method used 
for the determination of the band structure of the super- 
lattice. The following section discusses the assumptions 
made in the ultimate efficiency calculations and presents 
the method used for estimating the ultimate efficiency of 
a system with an intermediate band in the energy gap. 
Our numerical calculations are presented and discussed in 
a separate section. The paper is summed up in the Con- 
clusions and supplemented by an Appendix that presents 
Fourier coefficients of those material parameters which 
are periodic functions of the position vector (the effec- 
tive mass and the effective potential) for different rod 
geometries. 



II. SUPERLATTICE GEOMETRY AND BAND 
STRUCTURE 



Let us consider a two-dimensional superlattice formed 
by a system of GaAs rods disposed in sites of a square 
lattice and embedded in an AlGaAs matrix (see Fig. [1]). 
The rod axes are oriented along the (001) direction (the 
z direction) of the AlGaAs crystal lattice. The rods rep- 
resent potential wells for both electrons and holes propa- 
gating in the ix^y) plane perpendicular to the rod axes. 



a) 



b) 




c) 



d) 




Figure 1: Structure of two-dimensional semiconductor super- 
lattice, cross section perpendicular to rod axes (a nine-cell 
fragment of the infinite superlattice is shown). Black and 
white regions represent rod (well) and matrix (barrier) ma- 
terials, respectively. Unit cells are delimited by dashed lines. 
Centers of rods of each geometry coincide with centers of unit 
cells. In each depicted superlattice the filling fraction value 
is / = 0.25. Rods of cross section in the form of (a) an equi- 
lateral triangle, (b) a square and (c) a regular hexagon are 
oriented so that symmetry axes of their cross sections coin- 
cide with symmetry axes of the unit cell. 



The depth of the wells is determined by the concentra- 
tion of Al in the matrix material. Increasing Al concen- 
tration causes the effective potential felt by conduction 
band electrons or valence band holes to rise or decrease, 
respectively, with consequent deepening of the wells. 

Two-dimensional lattices formed by rods of different 
geometry are schematically represented in Fig. [TJ In this 
paper we consider rods of cross section in the form of an 
equilateral triangle (a), a square (b), a regular hexagon 
(c) or a circle (d). Rods are situated centrally in square 
unit cells of the superlattice and oriented in a manner 
which conserves as much as possible the symmetry of 
the lattice in the x and y directions. The filling fraction 
/, or the area ratio of the rod cross section to the unit 
cell, is a measure of the rod size. The maximum fill- 
ing fraction value corresponds to the situation in which 
the rods reach the unit cell borders. For rods of triangu- 
lar, hexagonal and circular cross-sectional geometry these 
maximum filling fraction values are 3\/3/16, 3\/3/8 and 
7r/4, respectively. 

The structure of energy minibands in the conduction 
band and in the valence band is determined here in the 
effective mass approximation, justified by the occurrence 
of a direct gap at point F of the atomic band structure 
for low concentrations of Al in both the rods and the 
matrix. To meet the conditions of applicability of the 
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effective mass approximation the Al concentration in the 
matrix material is only allowed to range from to 0.35. 

The following assumptions are made regarding elec- 
trons and holes: 

• Interactions between the conduction band and the 
valence band are negligible, allowing independent 
determination of the miniband structure for elec- 
trons and holes. 

• Only the interactions between heavy and light hole 
bands in the valence band are taken into account. 

• The z component of the wave vector is zero (elec- 
trons and holes propagate in the periodicity plane, 
(x,?/)), which implies a spin degeneracy of heavy 
and Hght hole bands. 

Electronic states are described by the Ben Daniel Duke 
equation^^: 
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with effective mass TO*(r) (isotropic in a homogeneous 
medium) and effective potential Ec{r) determining the 
position of the conduction band bottom; ^'e denotes the 
envelope wave function of a conduction band electron. 
The constant a = IQ-^^h'^ / {2mee) « 3.80998 (me and 
e denote the free electron mass and charge, respectively) 
allows to express the energy E in electronvolts (eV) and 
the coordinates x, y, z in angstroms (A). These mate- 
rial parameters vary in space with the periodicity of the 
superlattice: 



TO*(r + R) 
Sc(r + R) 



= TO*(r), 



(2) 



where R is a superlattice vector and r — (x, y) is the 
position vector. 

The Schrodinger equation of the envelope function of 
light and heavy hole states near the top of the valence 
band reads^^: 
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where ^'i;i(r) and "^hhit") are envelope functions for light 
and heavy holes, respectively, and operators P, Q and R 
have the form: 
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Luttinger parameters 71, 72, 73, describing, respec- 
tively, the effective masses 1/(71 -I- 72) and 1/(71 — 72) 
of light and heavy holes near the point F of the atomic 
lattice, are, like the position of the valence band top Ey , 
periodic in the superlattice structure. 



7/3(1- + R) = 7;3(r), 
Ev{y + K) = Evir), 



(5) 



where label /? denotes: 1,2,3. 

The effective mass of carriers and the effective poten- 
tial in which they move depend on the space-variable 
chemical composition of AlGaAs (alloy between GaAs 
and AlAs). Rods, which represent potential wells, are 
made of GaAs; the matrix material, representing a bar- 
rier for electrons and holes, is AlGaAs. The following 
empirical formulae, obtained by linear extrapolation of 
values of material parameters in GaAs and AlAs, allow 
to estimate their values in the matrix material^li^: 

Ec = 0.944d, 

Ev = 1.519-l-0.75d, 

m* = 0.067 -I- 0.083(i, 

71 = 6.85 - 3.40d, 

72 = 2.10-1.42d, 

73 = 2.90 - 1.61d, (6) 

where d is the concentration of Al in gaUium arsenide. 

The expansion in the plane- wave basis of the envelope 
function for holes and electrons and the Fourier expan- 
sion of the material parameters can be formally written 
as: 
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where the labels /3i, (32, refer to e, hh, Ih; C, V; 1, 2, 3, 
respectively; (j)'^^ have the sense of Fourier coefficients of 
the periodic factor of the envelope function (which has 
the same periodicity as the superlattice). Fourier coeffi- 
cients of the material parameters can be found analyti- 
cally from the formula: 
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for each of the four rod geometries under consideration. 
Symbols f{r) and F*^ denote, respectively, the periodic 
material parameter (effective mass or effective potential) 
and the corresponding Fourier coefficient for a plane wave 
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of wave vector equal to vector G of the reciprocal super- 
lattice. Parameter S represents the area of a superlattice 
cell. For the explicit form of the Fourier coefficients ([7]) 
for rods of triangular, square, hexagonal or circular cross 
section, see the Appendix. 

By substitution of expansions ^ and ^ in the 
Schrodinger equations |[T| and ^ for electrons and holes, 
respectively, we obtain systems of equations in the form 
of an eigenvalue problem, the solution of which yields 
Fourier coefficients of the periodic factor of the envelope 
function and the carrier energy. For electronic states we 
get: 



J2[a{G + k)- (G' + k) 
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In the case of interacting light and heavy holes the eigen- 
value problem reads: 
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where Ai,A2 and -Bi, B2 are expressed as follows: 



Ai = 
A2 = 
Bi = 

B2 = 



-a-ff-^' (G + k)- (G' + fc) - E^-^', 
-ajf-^' {G + k)- (G' + fc) , 
-aV37f [(G, + fc,)(G; + fc,) 

- {Gy+ky){G'y+ky)] , 

aVSj^-""' [{G^ + k^){G'y + ky) 

+ (G',+k^)iGy+ky)] 



(11) 



Figure [2] presents the band structure of a superlattice 
consisting of cylindrical rods (material A) embedded in 
a matrix (material B) with Al concentration d — 0.35, 
for a filling fraction / — 0.5 and a superlattice constant 
a — 8OA. The electron energy is referred to the energy 
EcA = of the conduction band bottom in the rod mate- 
rial; EcB is the position of the conduction band bottom 
in the matrix material; Eva and Eys denote the posi- 
tion of the valence band top in the rod and matrix mate- 
rials, respectively. The dispersion relation for minibands 
is plotted along the high-symmetry path T — X — M — T 
shown in the inset. 

In the (x, y) plane the mass of light and heavy holes 
is greater than that of electrons, which implies valence 
minibands narrower than conduction minibands. As the 
miniband structure is generated by both the light and 
heavy hole systems, valence minigaps are also narrower 
than conduction minigaps. Interactions between the sys- 
tems of light and heavy holes can cause additional narrow 
minigaps to open and minibands to split as a consequence 
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Figure 2: Energy spectrum of semiconductor superlattice 
formed by cylindrical rods of GaAs (A) embedded in a matrix 
of AldGai_dAs (d = 0.35) (B), with filling fraction / = 0.5 
and superlattice constant a = 80A. The dispersion relation 
is plotted along the F — X — M — F path shown in the inset. 
The horizontal lines indicate the position of the conduction 
band bottom in the rod and matrix materials {Eca and Ecb, 
respectively) and the position of the valence band top in the 
rod and matrix materials {Eva and Evb, respectively). The 
reference energy level is Eca = 0. The lowest miniband, de- 
tached from the conduction band, is regarded as an interme- 
diate level that opens an extra channel for carrier transitions 
between the valence band and the conduction band. The ul- 
timate efficiency of the solar cell is determined by the gap Eg 
between the conduction band and the valence band, and the 
distance Ej between the bottom of the intermediate miniband 
and the conduction band. 



of the lifting of the degeneracy of the two systems. Nu- 
merous very narrow valence minigaps are seen to occur 
in the well {Evb < E < Eva), and valence minibands to 
overlap below the barrier height {E < Evb) (Fig.E]). To 
simplify the ultimate efficiency calculation we assume the 
valence band is uniform (i.e. the minibands are merged) 
in the superlattice structures considered, and the highest 
miniband marks the top of the valence band. 



In a wide range of model parameters only the first (low- 
est) conduction miniband is distinctly detached from the 
band. This isolated miniband lies in the potential well 
{EcA < E < Ecb)- Minibands above the barrier height 
{E > Ecb) are inclined to overlap and tend to form a 
uniform block. 

In a very rough approximation, the band structure of 
the superlattice can be regarded as composed of a va- 
lence band, or block, (VB) and a block of overlapping 
conduction minibands, with a wide gap Eq between the 
two blocks. The lowest, detached conduction miniband 
can be regarded as an intermediate band (IB) within the 
gap. The energy Ej indicated in Fig. [2] is the shift be- 
tween the bottom of the IB and the block of overlapping 
conduction minibands. The latter shall be henceforth 
referred to as the conduction band (CB). 

If the symmetry of the rod cross section is lower than 
that of a square, points X and X' in the Brillouin zone, as 
well as points M and M', are not equivalent (see Fig. [21 
inset). In this study rods of triangular and hexagonal 
cross section are oriented against the superlattice in the 
manner shown in Fig.[TJ Propagation of an electron wave 
along the T — X direction is easily seen not to be equiva- 
lent to its propagation along T — X'. On the other hand, 
propagation along T — M and F — M' is identical. To 
examine the effect of the rod symmetry on the superlat- 
tice spectrum we have compared the dispersion relations 
for structures with rods of circular and triangular cross 
section. Figure [3] shows the dispersion relation along the 
path T — X — M — T — X' — M in superlattices with rods 
of both geometries. The calculations were performed for 
a filling fraction / — 3^3/16, the value corresponding 
to the triangular rods touching the borders of the super- 
lattice cell. In the superlattice with rods of triangular 
cross section the dispersion relation near points X and 
X' is only seen to differ visibly in higher minibands. Ob- 
viously, no such differences are found in the spectrum of 
the superlattice with cyHndrical rods. However, the two 
spectra do not differ much. The only differences are the 
lowest conduction miniband (the IB) slightly shifted to 
higher energy and the CB slightly shifted to lower energy 
in the superlattice with rods of triangular cross section. 



III. ESTIMATION OF ULTIMATE EFFICIENCY 
OF SYSTEMS WITH IB 

In our estimation of ultimate efficiency rj we made 
the following assumptions regarding the radiation of the 
source and the solar cell. 

• The cell absorbs the amount of black-body radia- 
tion at a temperature corresponding to the maxi- 
mum of the solar spectrum (the value assumed in 
the calculations is Tg = 5760^). This assump- 
tion simplifies our calculation and could affect the 
results only if optical transitions in the solar cell 
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Figure 3: Energy spectrum of conduction minibands of su- 
perlattices with rods of circular or triangular cross section 
(circles and filled triangles, respectively). Plotted along the 
path r~X-M-T-X'-M shown in Fig. El the disper- 
sion relation shows non-equivalence of Brillouin zone points 
X and X' in the superlattice with rods of triangular cross 
section. The following parameter values are assumed in both 
cases: superlattice constant a — 80A, matrix Al concentration 
d = 0.35 and filling fraction ft = 3\/3/16. 



structure coincide with windows in the terrestrial 
sun spectrum. 

• The temperature of the solar cell is T = OK, which 
impHes no recombination processes. The absence 
of recombination processes allows us to discard the 
impact of geometry of the solar cell and its align- 
ment with respect to the incident solar radiation. 

The assumptions listed below describe the features of the 
intermediate band photoconverter which is able to reach 
the maximal utilization of photon energy in the absence 
of carrier multiplication processes. The computed ulti- 
mate efficiency results only from selective absorption in 
a given band structure of the photoconverter. 

• Electronic transitions always occur through energy 
gaps of the greatest width available to the given 
photon energy. This implies minimized losses in 
the thermalization process. 

• If the photon energy is high enough for an electron 
to cross the gap, the probability of the transition is 
100%. The selection rules for optical transitions^^ 
are not taken into account. 

• The energy of excited carriers is not utilized until 
the electron crosses the gap to reach the conduction 
band as a result of cascade transitions via interme- 
diate bands. Then each electron brings an initial 
energy equal to the width of the gap Eq between 
the VB and the CB (Fig-EJ. 
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• The flux of electrons excited to the conduction band 
via intermediate bands is equal to the lowest flux 
of photons absorbed in a single gap in the cascade: 
VB-IB- CB. 

• Each absorbed photon generates a single electron- 
hole pair. 

A further increase of efficiency can by achieved by hot 
carrier extraction (when photon energy E > Eg) or by 
carrier multiplication (inverse Auger process). Hence, 
the assumptions presented above give the least restric- 
tive limitation for efliciency in the absence of carrier mul- 
tiplication and hot carrier utilization. The definition of 
ultimate efliciency we used results mainly from the photo- 
converter band structure and does not take into account 
the external geometry of the system (the fight manage- 
ment and concentration), the mismatch to the value of 
load (the shift of the electrochemical potentials), etc. 

The power of radiation per unit area of the solar 
cell is described by Planck's distribution for black-body 
radiatioiii: 



(12) 



where fcs is the Boltzmann constant, h is the Planck 
constant and c is the velocity of light. Three ranges of 
photon energy E can be distinguished in the considered 
case with a single IB. These three ranges correspond to 
the following three types of carrier transitions (Fig. [2]): 

• between the IB and the CB: Ei < E < Eq - Ei, 

• from the VB to the IB: Eg - Ej < E < Eg, 

• from the VB directly to the CB: E > Eg- 

The flux of photons absorbed by a unit area of the cell 
per unit time is governed by the following general formula 
for each of the above transitions: 



I{Ei,E2) = 2TT{kBT,)yh^c^ / 



- E/keT,, (13) 



where Ei and E2 denote the above-specifled limits of the 
energy range for the given transition. The utilized power 
per unit area of the ce\^: 



Pout =Eg[I (Eg, 00) 

+ min {IIEg - Ei,Eg),I{Ei, Eg - Ej j) 



(14) 



has two components. One, I {Eg, 00), is related to direct 
carrier transitions between the VB and the CB. The 
other component is the minimum value of two streams, 
I{Eg — Ej,Eg) and I{Ei,Eg — Ej), from cascade-like 
transition via the IB. The ultimate efliciency is deflned 
as a ratio of utilized output power to incident power: 
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Figure 4: Ultimate efficiency versus bandgap width Eg for 
two solar cells, without (solid line) and with (dashed line) an 
IB shifted by Ei from the conduction band edge. The ef- 
ficiency of the system with IB was determined for optimum 
position of the IB in the bandgap (see the inset), correspond- 
ing to equal fluxes of photons absorbed in transitions from 
VB to IB and from IB to CB. Indicated in the plot, the 
maximum efficiency values for the systems without and with 
IB are 42% (for Eg = 1.1 eV) and 68% (for Eg = 1.9 eV), 
respectively. 



Figure |4] presents the ultimate efficiency plotted versus 
the energy gap width Eg for two solar cefis, one without 
and the other with an IB in the bandgap. The assumed 
position of the IB corresponds to the situation in which 
the fiux of photons generating transitions between the 
VB and the IB is equal to the flux of photons generat- 
ing transitions between the IB and the CB. This means 
both fluxes are fully used for carrier transitions between 
the VB and the CB. The inset in Fig. |4] shows the op- 
timum relative position of the IB in the gap between 
the VB and the CB plotted versus the gap width Eq- 
The superimposed plot for the two systems under consid- 
eration indicates the ultimate efliciency is substantially 
higher in the system with IB, in which case the max- 
imum ultimate efficiency, about 68%, corresponds to a 
gap width Eg = 1.9 eV (in the solar cell without IB the 
maximum efliciency is 42% and corresponds to Eg = 1.1 
eV). 

These properties suggest AlGaAs could be a good can- 
didate for applications in solar cells based on materials 
with an IB in the bandgap, as the bandgap width in 
the monolithic materials ranges from 1.52 eV in GaAs 
to 3.21 eV in AlAs. However, in the superlattice model 
under consideration the IB, being the lowest conduction 
miniband, lies in the potential wefi between the conduc- 
tion band bottom in the rod material (GaAs) and in the 
matrix material (AlGaAs). Consequently, the distance 
Ej between the IB and the block of overlapping con- 
duction minibands is limited and far from the optimum 
value; this represents a major obstacle in attaining max- 
imum ultimate efliciency. However, the results presented 
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in the next section indicate a sensible increase in ultimate 
efRciency with respect to its values in the monolithic ma- 
terials without IB (AlGaAs or GaAs) in a wide range of 
superlattice geometry parameters. 



IV. RESULTS 

We have determined the ultimate efficiency for struc- 
tures depicted in Fig. [T] and examined it versus the fol- 
lowing structural parameters: the fiUing fraction /, the 
superlattice constant a and the concentration d of alu- 
minum in the matrix material, the latter parameter de- 
termining the depth of the effective potential well. The 
ultimate efficiency has been plotted versus each of these 
parameters, with the others fixed. We realize that a com- 
prehensive study should take into account the simultane- 
ous interplay of all mentioned parameters but such inves- 
tigation would be very difficult in analysis and require an 
enormous amount of computational time. The values of 
fixed parameters: / = 0.5, a = 80A , d = 0.35, were 
chosen arbitrarily. We tried, however, to give some con- 
vincing justification for such choice. The filling fraction 
/ = 0.5 describes the system in which the volume of the 
matrix and the rod rnaterials are equal; for the superlat- 
tice constant a ~ &0A the band spectrum looks quite reg- 
ular (the minibands and the minigaps are of comparable 
width); increasing the depth of the wells and the width 
of the IB (which results in higher efficiency) , the range of 
Al concentration d in AlGaAs is limited to d = 0.35 for 
the envelope function approximation to be applicable. 

Figure [5] shows the ultimate efficiency 77 versus the fill- 
ing fraction / for superlattices with rods of triangular, 
square, hexagonal or circular cross section. We fixed the 
superlattice constant at a = 80A and the matrix Al con- 
centration at d = 0.35. In the range of very low filling 
fraction values the effective potential wells formed by the 
rods become too narrow to bind electronic states. Pushed 
out of the well, the conduction miniband merges into the 
continuum of higher minibands; as a result, the superlat- 
tice band structure resembles that of monolithic AlGaAs 
with a continuous conduction band. The gap between 
the conduction band and the lowest miniband closes at a 
fining fraction value of 0.05, approximately. Note the ef- 
ficiency of the superlattice with a minimum width of the 
gap between the IB and the VB is relatively high, about 
31%, against the 26.8% ultimate efficiency of a homoge- 
neous matrix material (AlGaAs). This rather substantial 
difference can be explained as follows. When the effective 
potential well is narrow the position of the IB varies the 
most; the energy of the IB increases relatively fast with 
decreasing filling fraction. However, as the IB is pushed 
out of the well, it grows wider, and at the moment of clos- 
ing of the gap that separates it from the VB its width 
is quite significant and its bottom Hes below the effec- 
tive potential of the matrix. In this limiting case Ej (cf. 
Fig. [2I) is substantially different from zero, while Eq is 
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Figure 5: Ultimate efRciency rj versus filling fraction /, as 
calculated for superlattice constant a = 80A and matrix Al 
concentration d = 0.35. Each line refers to a different rod ge- 
ometry, as indicated by labels of respective shape. The hori- 
zontal line indicates the efficiency of monolithic rod material, 
GaAs. 

close to the width of the gap in monolithic AlGaAs. This 
affects significantly the ultimate efficiency in the model 
considered. 

The ultimate efficiency increases rapidly as the filling 
fraction continues to grow. At / = 0.18 the efficiency 
of a superlattice-based solar-energy converter is equal to 
that of a solar cell made of monolithic GaAs (the hor- 
izontal line in Fig. [5]). This steep increase in efficiency 
is due to the widening of the potential wells formed by 
the rods, and the concurrent narrowing of the barriers 
formed by the matrix material. As the wells grow wider, 
the first conduction miniband and the first valence mini- 
band move towards the well bottom, which results in a 
reduced Eg- Consequently, transitions omitting the IB 
become more intense, resulting in an increased efficiency. 
However, narrowing miniband spacing implies reduced 
distance between the IB and the CB. This competitive 
effect reduces the rate of transitions via IB. The ulti- 
mate efficiency reaches maximum at filling fraction 0.42, 
approximately. Further widening of the wells, or thick- 
ening of rods at the cost of the matrix material, leads to 
a decrease in efficiency. 

Note the filling fraction cannot reach 1 in structures 
with rods of triangular, hexagonal or circular cross sec- 
tion. Its Hmit values, 3-\/3/8 and 7r/4, for rods of hexag- 
onal and circular cross section, respectively, correspond 
to the situation in which the rods touch at the borders 
of superlattice unit cells. In the case of rods of trian- 
gular cross section the Hmit value / = 3-\/3/16 assumed 
in the calculations corresponds to rods touching the unit 
cell borders rather than one another, and is motivated 
by technical reasons and otherwise very complex calcu- 
lations. 

The minibands grow in width as the fiUing fraction 
continues to increase. In the superlattice with square 
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cross-section rods the gap between the IB and the CB 
closes at / = 0.81. The bottom of the continuous CB 
formed in this way is very close to the conduction band 
bottom in GaAs. However, the ultimate efEciency values 
of monolithic GaAs and a zero-gap superlattice differ by 
1%, approximately. This quite significant difference is 
due to the fact that as long as the gap remains open 
the conduction band bottom refers to the block of higher 
minibands. It is from this level that the energy of excited 
electrons is utilized in our model. When the IB merges 
into the CB the energy of utiHzed carriers is reduced by 
the width of the IB. 

As long as the filling fraction value is fixed, the ge- 
ometry of the rods proves of little importance for the 
efficiency. The efficiency is seen to increase as the shape 
of the rods becomes more compact and regular: the su- 
perlattice with cylindrical rods has a slightly higher effi- 
ciency than the superlattices with rods of triangular or 
square cross section, the approximate maximum differ- 
ences being 1.2% and 0.4%, respectively. 

Figure [6] shows the ultimate efficiency plotted versus 
the concentration d of aluminum in the matrix material 
for superlattices with rods of square, hexagonal and cir- 
cular cross section. The Al concentration in the matrix 
material determines the depth of the effective potential 
wells felt by carriers. The calculations were performed 
for the filling fraction fixed at / = 0.5 and superlattice 
constant a = 80 A. For d values around 0.14 the wells 
are too shallow to effectively bind electrons or holes, and 
the spectrum of the superlattice resembles that of mono- 
lithic GaAs. The intermediate band in these conditions 
is relatively wide, its bottom very close to that of the 
conduction band in monolithic GaAs, and the absolute 
gap between IB and the VB extremely narrow. The effi- 
ciency exceeds by 1% that of monolithic GaAs, which can 
be explained as in the case of the filling fraction depen- 
dence in the range of high filling fraction values. In the 
atomic structure of AlGaAs a direct bandgap only occurs 
at low Al concentrations. For the eflFective mass approxi- 
mation to be applicable we have confined the range of this 
parameter to d < 0.35. In the range 0.14 < d < 0.35 the 
ultimate efficiency increases almost linearly from 0.405 
to 0.425. This growth is due to the deepening of po- 
tential wells with increasing concentration of Al in the 
matrix. Deeper wells imply a greater distance between 
the IB and the VB formed by strongly overlapping hole 
minibands, which tend to lie above the matrix effective 
potential level. The efficiency increases with growing Ei 
(cf. Fig. [2]), due to an increased rate of transitions via IB. 
Carrier excitation in these cascade transitions is mainly 
hindered by a low fiux of absorbed photons able to induce 
transitions between the IB and the CB. 

Plotted in Fig. [71 the dependence on the superlattice 
constant a has been obtained for the filling fraction fixed 
at / = 0.5 and matrix Al concentration d = 0.35. The ef- 
ficiency is easily seen not to vary widely in the considered 
lattice constant range, 50A < a < 150A; the maximum 
difference is 1.2%, approximately. 
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Figure 6: Ultimate efficiency r\ versus Al concentration d in 
the matrix material (d determines the depth of potential wells 
formed by the rods with respect to the barriers formed by the 
matrix material). Calculations were performed for superlat- 
tice constant a — SOA and filling fraction / = 0.5. Plots 
obtained for superlattices with rods of square, hexagonal and 
circular cross section are labeled with respective geometric 
figures. 
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Figure 7: Ultimate efficiency rj versus superlattice constant 
a, as calculated for a filling fraction / — 0.5 and a matrix 
Al concentration d = 0.35. Plots obtained for superlattices 
with rods of square, hexagonal and circular cross section are 
labeled with respective geometric figures. 



Increasing the superlattice constant a at fixed filling 
fraction implies widening of both wells (associated with 
rods) and barriers (associated with the matrix). For 
small values of a the IB is weakly bound by narrow po- 
tential wells, in spite of the strong interactions between 
them. As a result, the IB lies beyond the well and is 
relatively wide. The distance Ej between the bottoms 
of the IB and the CB is relatively large; the decrease in 
efficiency is mainly due to the widening of the gap Eq 
between the VB and the CB. For large superlattice con- 
stant values the wells are wide and well isolated; in this 
situation both the miniband spacing and the miniband 
width are small; a number of non-overlapping minibands 
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can occur in potential wells. Our model assumes a sin- 
gle intermediate band. Omitting higher minibands that 
do not overlap can lead to an underestimation of the ef- 
ficiency. However, we believe that taking into account 
more minibands would not reverse the downward ten- 
dency of the efficiency with increasing superlattice con- 
stant. 

Note in the case of rods of circular and hexagonal cross 
section the efficiency reaches maximum at a superlattice 

constant of about 80^, while with square cross-section 
rods the optimum superlattice constant value is 95A. 



V. CONCLUSIONS 

We have estimated the ultimate efficiency of solar cells 
based on a two-dimensional superlattice of GaAs rods 
disposed in sites of a square lattice and embedded in an 
AlGaAs matrix. The efficiency of such superlattice-based 
solar-energy converters is found to exceed by a few per- 
cent that of monolithic GaAs in a wide range of param- 
eters. The key role in this gain in efficiency is played by 
the lowest conduction miniband, which, detached from 
the block of overlapping conduction minibands, acts as 
an intermediate band that opens an extra channel for car- 
rier transitions between the valence band and the conduc- 
tion band. The position of this intermediate miniband is 
determined by the distance Ej between its bottom and 
the bottom of the block of overlapping conduction mini- 
bands. Another parameter of vital importance for the ef- 
ficiency of solar-energy conversion is the distance between 
the top of the highest valence miniband and the bottom 
of the block of overlapping conduction minibands. This 
distance, or bandgap Eg, determines the energy of uti- 
lized carriers. 

The ultimate efficiency of a solar cefi with intermediate 
band is found to reach a maximum value of 68% for Eq = 
1.95 eV and Ei/Eq = 0.37. However, this maximum 



efficiency is unavailable in the material considered in this 
study, mainly due to the fimited distance Ej between the 
intermediate band and the conduction band from which 
it is detached. 

The filling fraction / proves to be the most important 
factor affecting the efficiency of the studied systems. For 
/ = 0.45 the ultimate efficiency exceeds that of mono- 
lithic GaAs by 3%, approximately. Of lesser impact is the 
superlattice constant, here allowed to vary in the range 
from 50 A to 150 A (at fixed /). Maximum efficiency is 
attained at a superlattice constant of about 90A. Inter- 
estingly, the geometry of rods is of negligible effect on the 
efficiency (as long as / is kept constant). The efficiency 
proves sfightly higher in superlattices with rods of more 
regular shape. 

Appendix A 



Defined as the area ratio of the rod cross section to 
the superlattice cell, the fifiing fraction {ft, fs, fh and 
fc, respectively) of superlattices with rods of triangular, 
square, hexagonal or circular cross section is expressed 
as follows: 

ft = ^{i)\ = (Al) 

where I, b and h denote the side length of the triangular, 
square and hexagonal cross section, respectively, and r is 
the cross-section radius of the cylindrical rods. 

Fourier coefficients F'^ (for G = 2n/2{m,n)) of a pe- 
riodic function /(r) (taking on the respective constant 
values /a and fs in the rods and the matrix) in a two- 
dimensional square lattice of rods of cross section in the 
form of an equilateral triangle (F^), a square (F^), a 
regular hexagon (F^) or a circle (F^) are given by the 
following formulae: 
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= < 



Fb- = < 



fB + ifA-fs)^ for m = n-0, 

if A - fB)^^e-'^-^'/^'' \^^^.nlVZ/a _ ^3^^^) _ 

for m = 0, n 7^ 0, 

if A - .fe) 2^.(^.l3,„„.) e-'=^""'/^a ^ g,:.„iV3/a (_73™cos(TO7r//a) 

+3m sin(m7rZ / a))] 
for m 7^ 0, n 7^ 0, 



/b + (/a - /b)^ 
(/>i-/s)^^sin(m7r6/a) 
(/a - /s);r^sin(n7r6/a) 



for TO = n = 0, 
for TO 7^ 0, n = 0, 
for TO = 0, n 7^ 0, 



F^ = 



(/a - fB):;^^^ [sin(TO7r6/a) sin(n7r6/a)] for m 7^ 0, n 7^ 0, 
/b + if A - /b)^#^ for to - n - 0, 

if A — fB) ^^^2^2 ['^ ~ a cos(\/3TO7r/i/a) + h^/Smn sini^/SmTrh/ ajj 
for m 7^ 0, n = 0, 

if A ~ f b) ^■n^rn'^n-ni) W^^ (cos(2n7r /i/a) - cos(\/3m7r/i/a) cosimrh/a)) 
+3m sin(\/3TO7r /i/a) sin(n7r/i/a)] 
for m 7^ 0, n 7^ 0, 

/b + if A - /B)7r^ for TO = n = 0, 

if A - /b) ^J2_^„2 Ji(27rVm2 + ?i2 ^/q) for to 7^, 71 7^0, 



(A3) 



where Ji(r) is a Bessel functfon of the first kind. 
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